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ABSTRACT. We consider gradient fields (4> x : x 6 Z d ) whose law takes the Gibbs-Boltzmann 
form Z -1 exp{— £\ B > V(cf>y — 4>x)} where the sum runs over nearest neighbors. We assume 
that V admits the representation 

V(n) = - log j £>(d«) exp [-|kt? 2 ] 

where q is a positive measure with compact support in (0, oo). Hence V is symmetric and non- 
convex in general. While for strictly convex V's the translation-invariant, ergodic gradient Gibbs 
measures are completely characterized by their tilt, a non-convex potential as above may lead to 
several ergodic gradient Gibbs measures with zero tilt. Still, every ergodic, zero-tilt gradient Gibbs 
measure for the potential V from above scales to a Gaussian free field. 



1. Introduction 

Gradient fields belong to a class of models that arise in equilibrium statistical mechanics, e.g., 
as approximations of critical systems and as effective interface models. Although their definition 
is rather simple, and in fact quite a lot is known (see the reviews by Funaki |fI51 , Velenik |[28l 
or Sheffield [25]), much is still there to be learned. In this note we study gradient fields on a 
lattice. Here the field is a collection of real- valued random variables = : x € Z } and the 
distribution of (j) on M z is given by the formal expression 

- exp{- ]T V(<j>y - J] d<fe, (1.1) 

(x,y) x&, d 

where d(f> x is the Lebesgue measure, (x, y) refers to an unordered nearest-neighbor pair on Z d , 
and V is an even, measurable function — called the potential — which is bounded from below 
and grows superlinearly at ±oo. 

Of course, to define the measure (11.11 ) precisely we have to restrict the above expression to a 
finite subset of Z rf and fix the </>'s on its boundary; Z is then the normalizing constant. Another 
way to regularize the expression (11.11 ) is to consider directly measures on all of R z whose con- 
ditional probabilities in finite sets take the above form. In our context this standard definition is 
hampered by the fact that, due to the unbounded nature of the fields <p x , no such infinite- volume 
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measures may exist at all. However, if one restricts attention to (the cr-algebra generated by) the 
gradient variables, 

Vxy ■= 4>y ~ 4>x, \x-y\ = l, (1.2) 

then infinite volume measures exist under the above conditions on V. Since the measure depends 
only on gradients, we refer to such measures as gradient Gibbs measures (GGM), in accord with 
Funaki El and Sheffield 1251 . 

Throughout we will focus on translation-invariant GGMs. An important characteristic is their 
tilt. For a translation-invariant GGM \i there exists a unique tilt vector t € M. d such that 

EM = t-b (1.3) 

for every edge b of 7L d — which we regard as a vector in this formula. Of course, this definition 
is really meaningful only for the GGMs that are ergodic — i.e., trivial on the cr-algebra of events 
invariant under all lattice translations. Indeed, in the ergodic case t represents the average incline 
of typical configurations. 

For the case of quadratic V — the massless free field — the measure (11.11 ) is Gaussian and 
so many desired characteristics are amenable to explicit computations. The challenge for math- 
ematicians has been to develop an equivalent level of understanding for non-quadratic Vs. A 
good amount of progress in this direction has been made in the last 10 to 15 years: Brydges 
and Yau (and also earlier works, e.g., by Gawedzki and Kupiainen lfl2l and Magnen and 
Seneor [19 ]) studied the effect of analytic perturbations of the quadratic potentials and concluded 
that the large scale behavior is that of the massless free field. Naddaf and Spencer ll23l proved 
the same non-perturbatively for strictly convex potentials V and GGMs with zero tilt. The cor- 
responding extension to non-zero tilt was obtained by Giacomin, Olla, and Spohn lfl3l . For the 
same class of potentials, Funaki and Spohn [16 ] proved a bijection between the ergodic GGMs 
and their tilt. Sheffield IT251 characterized translation-invariant GGMs by means of a Gibbs varia- 
tional principle and extended Funaki-Spohn's results to fields taking only a discrete set of values. 
We refer to the reviews by Funaki lfl5l . Velenik [28 ] and Sheffield [25] for further results and 
references. 

As a unifying feature, all the (non-perturbative) results mentioned are based on the strict con- 
vexity of the potential V — be it for the avail of the Brascamp-Lieb inequality lT23l [T3l [T6l . 
Helffer-Sjostrand random walk representation |[23l[T3l . coupling to Langevin dynamics |[T6l . and 
the cluster swapping-algorithm lf25l . One would naturally like to have a non-perturbative ap- 
proach that works even for non-convex potentials. With this motivation, Biskup and Kotecky JH 
recently studied the GGMs for non-convex V that are a log-mixture of centered Gaussians, 

V( V ):=-logj Q (dK)e-^\ (1.4) 

where q is a positive measure with compact support in (0, oo). For the simplest non-trivial case, 

q := pd Kl + (1 - p)5 K2 (1.5) 

with K\ S> k-2 > 0, it was shown that, in d = 2, there is a value p € (0, 1) at which one can 
construct two distinct, translation-invariant, gradient Gibbs measures of zero tilt. 

The relevant conclusion from Q for the general theory is that the one-to-one correspondence 
between ergodic GGMs and their tilt breaks down once V is sufficiently non-convex. The next 
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natural question is to understand what happens to the scaling limit. The purpose of this note is 
to show that, regardless of phase transition, for potentials of the form (11.41 ). every translation- 
invariant, ergodic GGM with zero tilt scales to a Gaussian free field. 

The proof is based on the fact — utilized already in [ 5 ] — that (11.41 ) allows us to represent 
every GGM as a mixture over Gaussian gradient measures with a random coupling constant n xy 
for each edge (x,y). Its covariance is simply the inverse of the operator 

(CJ)(x):= K xy [f(y)-f(x)} (1.6) 

y: |y-a:|=l 

where we take, once for all, K xy = K yx . The fluctuations in the Gaussian measure can be analyzed 
by invoking a random walk representation; C K is the generator of a random walk with symmetric 
random jump rates known, equivalently, under the name random conductance model. The name 
arises naturally from the electrostatic interpretation of this problem; cf. Doyle and Snell (9[, 
in which one views 7L d as a resistor network with conductance n xy — or resistivity 1 / K xy — 
assigned to an edge (x, y). As it turns out (see Lemma [3T2T >. if the initial GGM is ergodic then so 
is the law of the conductances. This makes homogenization a possible tool. 

Much work has been done in the past two decades on random walks with random conductances. 
For our purposes it suffices to invoke two known results: Kipnis and Varadhan's [ 18 ] invariance 
principle (i.e., scaling of the random walk to Brownian motion) and Delmotte and Deuschel's iflOl 
annealed derivative heat-kernel bounds. (Note that in the Helffer-Sjostrand random walk repre- 
sentation, as used in l23l [l3l . one also has to study a random walk in a random environment. 
However, this random environment fluctuates in time, while in our case it is static.) This handles 
the fluctuations of the field; an important technical issue is thus the control of the mean. This is 
where the zero tilt restriction comes in (see Lemma [3~4l Corollary 15. 81 and discussion in Sect. [6]). 

Note: While this paper has been in the process of revision, further developments occurred in 
the study of gradient models with non-convex potentials. Cotar, Deuschel and Muller (H have 
shown that for non-convex perturbations of potentials V where the size of non-convex region 
is small compared to typical fluctuations of the field, the conclusions are as in the convex case. 
(Their precise condition is a bound on the L -norm of the negative part of second derivative.) 
This is a high-temperature result; work in progress by Adams, Kotecky and Muller [ 1] addresses 
the low-temperaure case when non-convexities are allowed only sufficiently far away from the 
absolute minimum of V. Our contribution remains valuable despite these advances as it appeals 
to all temperatures, including those at which a phase transition occurs. 

The plan of this paper is as follows: In Sect. |2] we define precisely the concept of the gradient 
Gibbs measure and state our main theorem. In Sect. [3] we introduce the extended gradient Gibbs 
measures and characterize their conditional marginals. This will naturally lead to the aforemen- 
tioned connections with random walks in reversible random environments. To keep the main line 
of the argument intact, we first finish proving our main result in Sect. @] and only then expound 
on the random- walk connections in Sect. [5] Sect. [6] is devoted to discussion of the limitations to 
zero tilt and some open questions concerning gradient Gibbs measures. 
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2.1 Gradient Gibbs measures. 

As mentioned above, infinite volume measures on the field variables (4> x ) may not always exist 
— particularly, in sufficiently low dimensions. To make our statements uniform in dimension, 
we will focus attention on the gradient variables. However, not even that will be entirely straight- 
forward because the gradient variables satisfy a host of "hard-core" constraints which, in a sense, 
encapsulate most of the interaction. Since 77 is gradient, one has 

Vxi,X2 + Vx 2 ,X 3 + Vx3,X4, + Vx4,,Xi = 0, (2.1) 

whenever (x±, ■ ■ ■ , X4) are the vertices of a cycle in Z d of length four. We will often write % for 
the positively-oriented edge b in Z d . We will only work with positively oriented edges and will 
use B(A) to denote the set of such edges with both endpoints in the set A C Z d . 

The constraints (12.11 ) are implemented at the level of the a priori measure which is defined as 

follows: Fix a configuration 77 G R B ( zd ) that obeys dZlTl and, for A C Z d finite, let ^»( A ) c) be the 
Lebesgue measure on the linear subspace of configurations (rj' b ) such that r]' b = rjb for all b G" B(A) 
and that rj obeys the constraints (12- 1 b - Note that, if is a configuration such that rj xy = 4> y — <ft x 

for every nearest-neighbor pair (x,y), then za* (A)C is to within a normalization constant the 
projection to gradient variables of the Lebesgue measure on {4> x : x G A} subject to the boundary 
condition <p. 

Next we will give a precise definition of the notion of gradient Gibbs measure. For a finite 
A C Z d , consider the specification which is a measure in the first coordinate and a function 
of the boundary condition in the second coordinate, that is defined by 

7A(dTfe(A)l%(A)0 := ) ex p{ ~ Yl ^(W)} ^A B(A)C) ( dr ?B(A))- ^ 1 - T > 



(x,y) 
xeA, J/GAU9A 



Here Za(%(A) c ) is the normalizing constant. 

Definition 2.1 Let <f ffl(A ) = a({i] b : b G B(A)}). We say that a measure fx on M B ( zd ) is a 
gradient Gibbs measure if the regular conditional probability |<^b(A) c ) m an Y finite A C Z d 
satisfies 

^{-\^b(ay)(v) = 7a(-|%(a)0 ( 2 -3) 

for \i-a.e. 77. 

Most of this paper is restricted to translation invariant gradient Gibbs measures. To define the 
required notation, for each x G Z d let t x : M B ( zd ) — > R B ( zd ) be the "translation by x" which acts 
on configurations 77 by shifting the origin to position x, 

(T~x'r])yz • — 1]y+x,z+x • (2-4) 

We say that p is translation invariant if p o t x 1 = p for all x G 7L d , and that it is ergodic 
if p(A) G {0, 1} for every event A such that t~ 1 (A) = A for all x G 1 d . 
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2.2 Scaling limit. 

As is usual, we will interpret samples from gradient Gibbs measures as random linear functionals 
on an appropriate space of functions. Let C^°(M. d ) denote the set of all infinitely differentiable 
functions /: M d — > E with compact support. Given a configuration r\ = (%) of gradients 
satisfying the conditions (12.11) . we can find a configuration of the field (ft = {(ft x ) such that (11.21 ) 
holds for every nearest-neighbor pair of sites. The configuration (ft is determined uniquely once 
we fix the value at one site, e.g., 4>q. For any function / € C^°(R. d ) we introduce the random 
linear functional 

<f>(f) := Jdxf(x)4> [x] (2.5) 

which, under the condition 

fdxf(x)=Q, (2.6) 

does not depend on the choice of the special value (fto- 

The functional (ft(f) can be naturally extended to a somewhat larger space defined as follows. 
Let A denote the Laplace differential operator and consider the set 

H :={Ag: geC^(R d )}. (2.7) 

Note that each / G Hq automatically obeys (12.6b - The set Ho is endowed with a natural quadratic 
form / i-> (/, /) + (/, -A" 1 /) defined as 

(Ag,Ag) + (Ag,-A- 1 Ag)= [ dx {\Ag(x)\ 2 + \Vg(x)\ 2 ) . (2.8) 

We thus define the norm 

\\f\\n ■= [(/,/) + (/, -A" 1 /)] 1/2 (2-9) 
and let Ti be the completion of TLq in this norm. Note that H is the case of k = — V2 in the family 
of Sobolev spaces Ty fc ' 2 (R d ). The condition (/, -A" 1 /) < 00 is natural once we realize that 
this quantity will represent the variance of the limiting Gaussian field. 

The extension of <j) to Ti is implied by the following lemma: 

Lemma 2.2 Suppose g in (11.41 ) has support bounded away from zero and let [ibe a translation- 
invariant, ergodic zero-tilt gradient Gibbs measure for the potential V. Then there exists a con- 
stant c < 00 such that for each f £ TLq, 

WU)\\l*m <4f\\n- (2-10) 
In particular, (f) extends to a linear functional (ft: Ti — > R. 

Note that (12.101 ) means that the map / 1— > <ft{f) is continuous in L 2 -norm. If we want to avoid 
questions about accumulations of null sets, this permits us to work with only a countable number 
of /'s at each time. (In particular, we do not claim that / 1— > cft(f) is continuous in any pointwise 
sense.) This will not pose any problem because we will content ourselves with the (weaker) 
definition of a Gaussian free field based on the standard approach via Gaussian Hilbert spaces 
(cf. Sheffield [26, Sect. 2.4]): 
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Definition 2.3 We say that a family {ip(f) ■ f £ T~t} of random variables on a probability space 
(CI, & ' , P) is a Gaussian free field if the map f i— > tp(f) is linear a.s. and each ip(f) is Gaussian 
with mean zero and variance 

E(^(f) 2 ) = (f,-A- 1 f). (2.11) 

Our goal is to show that the family of random variables {4>(f) : / £ TL} has asymptotically, in 
the scaling limit, the law of a linear transformation of a Gaussian free field. To pass to this limit, 
we have to impose that the test functions are slowly varying, which we take to be on the scale 
e~ x . For e > and a function f:R d ^R, let 

f £ (x) := 6^ +1 )/(6x) (2.12) 
and note that the normalization ensures 

||/ £ |^ = (/ e ,/ e ) + (/ (-A)- 1 /e)=e 2 (/,/) + (/,(-A)- 1 /) < \\f\\ 2 n- (2-13) 

e<l 

Let cf) e denote the linear functional acting on / £ (M. d ) via 

Mf) ■= ttfe) = j dxf(x) (e^^Ul./ei). (2-14) 
The main theorem is the Gaussian scaling limit for </> e (/): 

Theorem 2.4 (Scaling to GFF) Suppose V is as in (11.41 ) with g compactly supported in (0, oo). 
Let /i be a gradient Gibbs measure for the potential V which we assume to be ergodic w.r.t. 
translations and to have zero tilt. Then for every f £ Tt, 

Hm£ M (eW)) = exp{i J dxf(x) (Q~V)0)}, (2.15) 
where Q~ 1 is the inverse of the operator, 

o/'-Ewa^/ <216) 

with (qij) denoting some positive semi-definite, non-degenerate, d x d matrix. In other words, 
the law of (j) t on the dual £■' of any finite-dimensional linear subspace 6 C 7i converges weakly 
to that of a Gaussian field with mean zero and covariance (— Q) _1 . 

Remarks 2.5 Here are some additional observations and remarks: 

(1) Since (— Q) is dominated by a multiple of (—A) from below, the integral in (12.151 ) — inter- 
preted as the quadratic form (/, Q~ l f) — is well defined for all / £ Tt. 

(2) Note that, in (12- 14b . the individual </>'s get scaled by /2+1 ) and not f ' 2 as one might 
expect from the conventional central limiting reasoning. This has to do with the fact that 
the variables (<p x ) are strongly correlated. These correlations are weaker for the gradients 
Vxy = <f>y — <j>x which adhere to the "usual" central limit scaling. In d = 1 and for gen- 
eral potentials V, the increments 775 are in fact i.i.d. and the scaling limit follows from the 
standard CLT 
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(3) In d > 1 the matrix (g^) is not necessarily a multiple of unity, since, in general, [i is not 
guaranteed to be invariant under reflections and rotations of 7L d . (Nevertheless, we expect 
that an ergodic translation invariant measure for the isotropic interaction (12.21) will inherit 
these symmetries.) To get convergence of 4> € to GFF in the sense of Definition 12.31 one thus 
has to scale the argument of <fi by the root of the corresponding eigenvalue of q in each of its 
principal directions. 

(4) The absence of strict convexity does not permit us to use the general argument of Funaki and 
Spohn |[T6l for the existence of an ergodic GGM with zero (or any other prescribed) tilt. To 
show that such GGMs do exist — and that our Theorem I2.4l is not vacuous — we note that, 
by Lemma 4.8 of Biskup and Kotecky (5), every weak limit of torus measures exhibits expo- 
nential concentration of the empirical tilt; one then just needs to pick any ergodic component. 
Note that this lemma applies only to zero tilt (cf. ||5] Remark 4.9]). 

(5) The restriction to zero tilt is actually a significant drawback of our analysis. The main reason 
is our inability to characterize the scaling limit of the so called corrector for the corresponding 
random walk problem. See Sect.[6]for more details. 

(6) In the example studied by Biskup and Kotecky Q, cf. (11.51 ). the two GGMs coexisting at the 
transitional value p t of p were proved to exhibit different characteristic fluctuations. It follows 
that the corresponding scaling limits will be distinguished by their stiffness coefficients q^j. 

Moreover, by Theorem 2.5 of 0, for k\ S> &2 the transition in d = 2 model with (11.51 ) 
lies on a self-dual line, i.e., 



The transition presumably stays on this line even as one slides the ratio K1/K2 towards one. 
However, it disappears before K\ / K2 hits one because for K\ ~ ki the potential V is convex 
and so there is only one GGM with zero tilt iPToll . At such point of disappearance physicists 
often expect non-trivial critical fluctuations. Notwithstanding, our results show this is not the 
case. 

(7) We avoid the context of the "stronger" definition of GFF as a random element in an appro- 
priate Banach space (cf. Gross [14] or Sheffield [26, Sect. 2.2]). This definition is appealing 
in d = 1, where the limiting functional / i— » tp(f) actually admits the integral representation 



with t i— » ipt denoting a continuous diffusion with generator Q, but in d > 1 the correspond- 
ing field becomes less and less regular with increasing dimension and the appeal is lost. 
However, this context would be in demand if one wishes to discuss the notion of tightness 
and convergence in law for the limit in Theorem l2.4l 

Both Lemma [2721 and Theorem 12.41 are proved in Sect. [4] 




(2.17) 




(2.18) 
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3.1 Coupling to random conductance model. 



The key idea underlying the representation dl.4| ) is that the auxiliary variable k in the expression 
for V may be elevated to a genuine degree of freedom associated with the corresponding edge. 
Specifically, given a gradient Gibbs measure p with potential (1 1 -4b . for each finite A C B(Z rf ) 
consider the measure p\ on R B ( Z ) x M A defined by 

p A (A x B) := I H g(d Kb ) fiiufj e v{m,)-farA > ( 3. 1} 
beA ^ beA ' 

where A C M B(zd) and BcK A are Borel sets. The representation (11.41) ensures that (/2a) is a 
consistent family of measures; by Kolmogorov's Extension Theorem these are projections from 
a unique measure p on configurations (rjb, kj,) G M B ( zd ) x M B ( zd ). The restriction of /2 on the 
?]'s gives us back p; we call /2 an extension of /x. The measure p is Gibbs for the Hamilton- 
ian Y^( x y ) \ K xy r Jxy> so we w ^ re f er to it as an extended gradient Gibbs measure (see Biskup and 
Kotecky Q for further facts on the extended GGMs). 

To ease the notation, whenever b is an edge between x and y, we may interchangeably write 
Kb and K xy for the same quantity. Furthermore, for the same reasons it will even be convenient to 
assume 

^xy — K-yx, \x y \ — 1, (3-2) 
and work with the k's as symmetric objects. 

We proceed by a series of lemmas that characterize the properties of ft. 

Lemma 3.1 Let \i be a gradient Gibbs measure for the potential V and let \x be its extension 
to M B ( zd ) x M B ( zd ). Consider the a-field S := a({r/ b : b G B(Z d )}). For fx-a.e. r\, the regular 
conditional distribution \<f>)(rj), regarded as a measure on the k's, takes the product form 

fi(dK\£)( V ) = (g) [e v ^)-s^e(dK 6 )]. (3.3) 



Proof. Recall that &wa) := a ({ r lb : b G B(A)}). The identity (13.11 ) implies that /2a coincides 
with jl on M B ( A ) x M B ( A ). But fj,\(— \&\) has the desired product form by definition and so the 
claim follows by standard approximation arguments. □ 

Lemma 3.2 Let p, be a gradient Gibbs measure and let p be its extension to IR B ( zd ) x M^^ d \ 
If p is translation-invariant and ergodic then so is p. 



Proof. The uniqueness of the extension of measures (13.11 ) implies that p is translation-invariant 
if p is translation invariant, and so it remains to prove that ergodicity is inherited as well. Let 
A C M B ( zd ) x R B ( zd ) be a translation invariant event, i.e., (rj, k) G A iff (r x rj, t x k) € A for all x. 
Our task is to show that p(A) G {0, 1}. 

First we invoke the ergodicity of p. Consider the function 

f( V ):=Ez(l A \g)(r]). (3.4) 
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Since A and fi are translation invariant, we have 

f(r xV ) = Efi(l A \£)(t xV ) = ^(l A or-V)W = f(v), £-a.s.. (3.5) 

But / is S -measurable and the restriction of ft to $ is fj,, which we assumed to be ergodic. 
Hence / is constant almost surely. Let c denote this constant. 

We will use a standard approximation argument to show that c G {0, 1}. Since A is an event 
from the product cr-algebra, there exists a sequence of events 

A n G a({r] X y, K xy : \x - y\ = 1, [ar| < n}) (3.6) 

such that 

fl(AAA n ) — > 0. (3.7) 



n— >oo 

The bound 



II U, - 1a || z i(A) < HAAA n ) (3.8) 
shows that then l An — > 1 A in L (ft). Since A is translation invariant, we have 1a = 1a l T -i( A )- 
Each indicator can be approximated by the indicator of the event A n ; a simple bound gives 

|| l An l T -i (An) - 1 A 1 t - Ha) |Ui W < 2f,(AAA n ). (3.9) 

For x with |s| > 2n + 1, the fact that p,{—\S) is a product measure (cf Lemma |3~TI ) implies 
that A n and r~ 1 {A n ) are independent. Hence, 

E^l An \-, {An) \S) = E fi (l An \*)Ep(l T ; HAn) \*). (3-10) 

Rolling the approximations backwards, we then conclude that the left-hand side converges to c 
in while the right-hand side converges to c 2 (note that all expectations are bounded). It 

follows that c = c 2 , i.e., c G {0, 1}. As 

fl(A) = Efi(f) = c, (3.11) 

the proof is finished. □ 

3.2 Random walk connections. 

Our next goal will be to characterize also the conditional measure given the re's. This will in 
turn require some facts from the theory of random walks with random conductances. We will 
frequently borrow facts from an associated potential theory which will be expounded in Sect. [5] 

Pick a configuration k = (Kb) with «{, G (0, oo), and recall the formula (11.61 ) for the genera- 
tor C K of the random walk among conductances k. We will focus on the action of C K on functions 
of both the environment k and the position x that satisfy the following shift-covariance property 

g(K,x + b) -g(K,x) = g(r x K,b), (3.12) 

with x G Ij d and b a coordinate unit vector in U. d , subject to the condition 

g(K,Q)=0. (3.13) 

This makes the function completely determined by its values at x being the neighbors of the 
origin. A function of this kind is said to be harmonic for the above random walk if 

-) = (3.14) 
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for (almost) every k. As it turns out, harmonic, shift-covariant functions are uniquely determined 
(a.s.) by their mean with respect to ergodic measures on the conductances: 

Lemma 3.3 Let u be a translation-invariant, ergodic probability measure on configurations 
K = ( Kb ) € M. M ^ such that u(e < K b < Ve) = I far some e > 0. Let g : M B ( z<i ) xZ^Rka 
measurable function which is 

(1) harmonic in the sense of (13.141) , v-a.s., 

(2) shift-covariant in the sense of H3.12\\3.1 J I ), v-a.s., 

(3) square integrable in the sense that E u \g(-, x)\ 2 < oo for all x with \x\ = 1. 
If E u (g(-,x)) = Ofor all x with \x\ = 1, then g(-, x) = a.s. for all x G Z d . 

We defer the proof, and further discussion of the consequences of shift-covariance and har- 
monicity, to Sect. [5] Returning to the gradient fields, we now characterize the conditional law 
given the k's: 

Lemma 3.4 Let [xbe a translation-invariant, ergodic gradient Gibbs measure with zero tilt and 
let (I be its extension to x R M( * zd \ Consider the a-field & := a({K b : b G B(Z d )}). 

For pi-a.e. k, the conditional law fi(— regarded as a measure on the set of configurations 

{((ftx) € M.^ : 4>o = 0} with the (j)'s defined from the rj's via (11.21 ). is Gaussian with mean zero 

Efi(<l> x \&)(K) =0, xeZ d , (3.15) 

and the covariance given by (—C^ 1 . Explicitly, for each f : Z d — > M. with finite support and 

= o, 

Var A ^/(x)^ A(k) = Y J f(.x)(-C' 1 f)(x). (3.16) 

^ X ' x 

Proof. The fact that the conditional measure is multivariate Gaussian law with covariance £~ x is 
checked by direct inspection of (I3.ll ). The only non-trivial task is to identify the mean. First we 
note that the loop conditions (12.11 ) ensure that there exists a function E such 

that 

u(k,0)=0 (3.17) 

and 

u(k, x + b) - u(k, x) = Efi(r] XtX+b \^)(K) (3.18) 

for all unit vectors b in the coordinate directions. We claim that u is harmonic in the sense of 
(l3~T4l) . Indeed, 

£ k u(k, x) = Efi ( ^2 K xy ((py - 4> x ) ^ (k) (3.19) 

y: \y-x\=l 

where we write, thanks to the loop conditions, rj xy = <p y — (j) x . Using that fi is Gibbs, we can now 
also condition on <r(4> y : y ^ x); the conditional measure fi{ x y is Gaussian with the explicit form 

fi{x}( d ^>x) = ^ eX p{ _ 0^ ^ K X y + <f> x ^ K xy<Py^d<p x , (3.20) 

y: \y-x\ = l y: \x-y\ = l 
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where Z is an appropriate normalization constant. It is easy to check that the mean of <\> x Y^ y . ij,_a;|=i K - 
under fi {x} is exactly Y, y - \ y - x \=i K x y 4> y , proving C k u(k, x) = 0. 
Next we observe that the translation invariance of jj, implies 

u(r x K,b) - u(t x k,0) = Ejx(r] 0ib \J?)(T x K) 

= Eji(T} X)X+b \&)(K) =u(k,x + b) -u(k,x) (3.21) 

and so u is shift-covariant as defined in (13. 12j - [3~T3l) . Finally, the definition of u and the fact that fx 
has zero tilt imply 

Efi{u(',x)) = EpSjfyp) = 0, |x| = 1. (3.22) 
As u obeys all conditions of Lemma [331 we have E^{(\) x \,^) = u(-,x) =0 /i-a.s.. □ 

Our reference to the random walk with generator C K is not limited to Lemma [331 we will also 
need to know some specific properties of this random walk. First, we will need to know that the 
position of the walk satisfies a Central Limit Theorem. Let X = (X t ) denote the continuous- 
time random walk with the generator C K and let P% denote the law of X subject to the initial 
condition P%(Xq = x) = 1. The following lemma goes back to Kipnis and Varadhan |[T8l . 

Lemma 3.5 (Annealed CLT) Let p,be a measure on IR B ( z<i ) which is translation-invariant and 
ergodic, and obeys fi{e < uj^ < Ve) = 1 f or some e > 0. Then there exists a positive semi- 
definite, non-degenerate, d x d matrix q such that, for every t > 0, the annealed distribution 
E^P®(eX te ~2 6 •) converges weakly to the law of the multivariate normal M(0, tq). 

The main result of |[T8l actually shows that the annealed law of the entire path t i— > eX te -2 
converges to that of (a linear transform of) Brownian motion. However the above is all what will 
be needed for the purposes of the present paper. 

Apart from a central-limit asymptotics, we will also need to have an estimate on the heat-kernel 
of the above random walk. The following lemma is a consequence of the main result of Delmotte 
and Deuschel ifTUl . 

Lemma 3.6 (Heat-kernel upper bound) Let fi be a law of the conductances satisfying the ellip- 
ticity condition p,(e < Kb < Ve) = 1 for some e > 0. Then there is a c\ < oo such that 

Ep\V i V j P%(X t = >)\<-^, l<i,j<d,t>0. (3.23) 
Here Vj is the discrete derivative in the i-th direction, Vj/(x) := f(x + 8j) — f(x). 
Proof By formula (1.5b) in flU Theorem 1.1], 

E^ViVjP^Xt = x)\ < c[ Y~L, (3.24) 

where P c '2*(0,x) is the probability of the continuous-time simple random walk to be at x at 
time c' 2 t. This probability is bounded by a constant times t~ /2 . □ 



4. Proof of main result 
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4.1 Regularity estimates. 

The goal of this section is to prove Theorem 12 .4 1 concerning the scaling limit of <p(f e ). We begin 
by proving L 2 -continuity of the random functional / i— > 4>(f) on H as stated in Lemma I2T21 For 
convenience of notation, whenever 1Z is an operator on L 2 (Z d ) we will extend it to an operator 
on L 2 (R d ) via the formula 

(f,Tif) := J dxdyf(x)f(y)n{lx],[y]), (4.1) 

where lZ(x, y) is the kernel of 1Z in the canonical basis in L 2 (Z rf ). 

Proof of Lemma |Z2] Let jl be the extended gradient Gibbs measure corresponding to fj,. Recall 
the notation L K for the generator of the random walk among conductances k = (/«&) and let £ 
denote the generator of the simple-random walk (i.e., the special case of C K when all kj, = 1). 
Pick / € {Ag: g G Co°(R d )}. Lemma 1341 and the fact that <j>(f) is hnear in the 77's imply 

U(f)\\hw = M(f-£« 1 f)) ( 4 - 2 ) 

where (/, C^ 1 /) is as defined above. By assumption on the support of g we know that > a 
/i-a.s. by which we immediately have the operator inequalities 

(-£«)> a{-C) and {-C K )~ X <cr\-C)- x . (4.3) 

Therefore, it suffices to bound the quadratic form associated with the (homogeneous) discrete 
Laplacian C by the quadratic form defining the 7Y-space: 

(/,(-£)-7) <c\\f\\ 2 n (4.4) 

for some constant c < 00 and all / in a dense subset of Ti. 

To this end we pick / G H in the Schwartz class and let f(k) := (27r) -d / 2 J f(x)e. lk ' x dx be 
its (L 2 -norm preserving) Fourier transform. A direct calculation now yields 



fc':3fgZ d J =1 
k-k'=2nl 



(4.5) 



where 

d 

(-C)(k) :=4^sin 2 (A ;j /2) (4.6) 

3=1 

is the generalized eigenvalue of the lattice Laplacian. Introducing —A(k) = \k\ 2 to denote the 
corresponding quantity for the continuum Laplacian, we invoke the Cauchy-Schwarz inequality 
for the sum over k' to get 

(M-C)- 1 /) < f dk\ Yl l/WU-^'HW) (4.7) 

k-k'=2n£ 
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where 

c(k):= E i -Am fj( 2si y y (4 . 8) 

k-k'=2-Kl 

is well defined on the set B := {k G [— 7T, ir] d : kj 7^ 0, j = 1, . . . , d} of full Lebesgue measure 
in [—it, ir] d . We claim that 

c := sup c(k) < 00. (4.9) 
fces 



We will show this by proving that the summand in (14.81 ) is bounded by a constant times the 
product rijd^jl + l) 2 - Indeed, for the A;' = A; term we use that the ratio A (k)/C{ k) is bounded 
throughout B, and same for the ratios k 1— > svn{kj/2)/kj. When At' 7^ fc, we set i to be the first 
index j where k'- 7^ kj and bound the 4sin(A^/2) 2 term by —C(k). Then we bound the ratio of 
A (A;) -terms by unity and the j-th term in the product by a constant times (|A^| 2 + l) -2 - (The 
sin {k'i/2) -term is not needed because \k[\ > 2tt.) 

The product PX/d^jl + ^) 2 * s summable over k! G k + (27rZ) d uniformly ink £ B and so 
(14.91 ) is proved. Bounding c(A;) by its supremum in (14.71 ). we can merge the sum and the integral 
to get ||/ 1 1 2 ^. The desired bound (1441 ) follows. □ 

Remark 4.1 The inclusion of L 2 -norm of / into \\f\\n is crucial for the bound ( I4.4I ). Indeed, on 
the basis of (14.51 ) it is not hard to construct functions for which the ratio (/, —C^ 1 /) / (/, — A -1 /) 
is arbitrarily large. This is caused by the fact that the spectrum of —A extends all the way to 
infinity while that of — C is bounded. 



The continuity established in Lemma [2721 allows us to work only with smooth and compactly- 
supported test functions. We will nevertheless need one more regularity bound before we can 
delve into the proof of our main result: 

Lemma 4.2 Let n be a translation-invariant law on the conductances subject to the ellipticity 
condition, fi(e < < Ve) = 1-for some e > 0. Then there exists c < 00 such that whenever 
f = Ag for some g G Cg°(R d ), 

E^f^f) < c\\Vg\\l A(su PP5 ) 2 ^_, (4.10) 

where X(A) is the set function on Borel subsets ofM. d defined by 

d . 

A (^) : = EE/ dzl {*+*eA} (4-11) 

^^Ud i = 1 



with Ri denoting the set of points in [0, l] d with vanishing i-th coordinate. 
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Proof. Translation invariance of \i and a simple integration by parts tells us 

E^Ag^Ag) = V / dz f dz' Ag{x + z)Ag(y + z')E^P% (X t = y) 



= J2 [ dz [ dz'd i g(x + z)d j g(y + z')V i V j E^(X t = y), (4.12) 

x,y&Z d i,j=l jRi R i 

where 8$ is the partial derivative with respect to the z-th coordinate. Restricting the integrations 
and sums so that the arguments x + z and y + z' are in the support of g, bounding the partial 
derivatives by [| V<?||oo> an d applying (13.231) . we obtain the desired bound. □ 

The consequence of Lemma l4~2l that will concern us is as follows: 
Corollary 4.3 For \i as in Lemma \4~2\ and any f £ {Ag: g £ Cg°(M d )}, 

/•CO 

lim sup / dt ^e- 2 (/ £ ,e te " 2£ «/ e ) = 0. (4.13) 

Proof. Pick / of the form / = Ag and note that f e = Ag( e \ where g( e \x) = e d/2+l g(xe). First 
we observe 

llV^Hoo = e^HV^IIoc. (4.14) 

Next we note that, since the support of g is the closure of a non-empty bounded open set, a simple 
covering argument tells us that 

e d A(supp 5 (e) ) d|suppg|, (4.15) 

where Isuppgj is the Lebesgue measure of suppg. As a consequence, there exists a con- 
stant C(g) < oo such that 

A(supp 5 (e) ) < C{g)e- d , < e < 1. (4.16) 
Plugging (14.141 ) and (14.161 ) into (14.101) . we get for e £ (0, 1), 

E M e- 2 (/ e ,e fc - 2 ^/e) < cIlVsH^) 2 ^. (4.17) 
The functions on the left (indexed by e) are uniformly integrable in t in all d > 1. □ 

4.2 Scaling limit. 

Having dispensed with regularity considerations, we can now proceed to establish the principal 
fact underlying the proof of Theorem l2.4t 

Proposition 4.4 Let fx be a translation-invariant, ergodic measure on k = (K b ) € R M ^ such 
that n(5 < Kb < Va) = 1 far some 5 > 0. Then there exists a positive semidefinite, non- 
degenerate, d x d matrix q = (g^,- ) such that 

lim(/ e , (-40- 1 /*) = (/,(-Q)-7) (4-18) 
40 
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in ^-probability for each f G C^°(M. d ) fl TC, where Q is defined from q by ( I2.16I ). 

Key to the proof is the following lemma: 
Lemma 4.5 For every t > and any f G C$°(R d ) n H 

&e(t) := e- 2 (/e,e' e " 2£K A) - (/,e«/) in L 2 {p). (4.19) 

Proof. Let (-Xt)t>o be the continuous-time random walk with the generator C K and let denote 
the law of the walk started from x. By Lemma [331 the annealed law of eX te -2 tends weakly to 
that of multivariate normal 

M t :=M{0,tq) (4.20) 

for some positive semidefinite, non-degenerate, d x d matrix q = (%•). As a consequence, if 
Q C Co°(M d ) is a family of functions that are uniformly equicontinuous and bounded, we have 

eJ sup E° K (g(eX te - 2 )) - Eg{M t ) *) — ► 0. (4.21) 

Now let Q be the generator of the Brownian motion with mean zero and covariance q; i.e., Q is 
the operator in (12- 16b - Then we have 

(/,e tQ /) = J dyf{y)E{f{y+N t )). (4.22) 

As C K is the generator of the random walk (X t ), we derive similarly 

e- 2 (/ e ,e fc_a£ -/ e ) = e d J Az x Az 2 ^ f(ex + e Zl )E*(f(ez 2 + eX te - 3 )). (4.23) 

[0,l] d x[0,l] d aeZ<i 

We thus need to show that the right-hand side of (14.231) tends to that of (14.22b - Note that if / is 
supported in [— M, M] d , the integral in (14.221) can be restricted to y G [— M, M] d and the sum 
over x in (14.231) to, say, \x\ < 2M/e (once e <g; 1). 

Substituting y = ex + ez with x G 7L d , \x\ < 2M/e, and z G [0, l] rf in (14.221 ) allows us to put 
both terms on the same footing. Subtracting (14.221 ) from (14.23b . taking expectation with respect 
to n and applying the Cauchy-Schwarz inequality we thus get 



E^itf < \\f\\ 2 e d J dz Yl E^\E x J{ez + eX tt -2)-Ef\ex + ez + Nt, 



2 . (4.24) 



[0,1]° 



x: \x\<2M/e 



Using translation invariance of \i, we may replace E*f(ez + eX te -2) by E®f(ex + ez + eX te -2) 
inside the expectation. For / G Co°(R d ), 

Q = {/(ex + ez + ■) : \x\ < 2M/e, z G [0, l] d } (4.25) 

is an equicontinuous family of uniformly bounded functions. Equation ( 14.211 ) then tells us that 
the right-hand side of (14.241 ) tends to zero as e [ 0, proving (14- 19b - □ 
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Proof of Proposition \4~4\ To extract d4~T8T > from d4~T9l , we note that for any / <G Cfi°{R d ) n H, 

{f^-C K )- l f) = / dt(/,e^/)- (4-26) 
J o 

Replacing / by f e and scaling t by e 2 , we find 

/■oo 

(/ e ,(->C re )- 1 /e)= / dt6- 2 (/ e ,e'^/ e ). (4.27) 

JO 

By ( 14. 191 ) the function under integration tends to (/, e*^/) in probability for each t; the mono- 
tonicity in t (and continuity of the limit) ensures that the convergence is actually uniform (in 
probability) on compact intervals. By Corollary 14.31 the integral can be truncated to a finite in- 
terval in L 1 -norm, and similarly for the integral of the limit which is finite as / is in the domain 
of Q- 1 . It follows that 

(/,,(-£„)-!/*)—► / dt(/,e«3/) = (4.28) 

in //-probability (and L l (//)). This is the desired conclusion (14.181) . □ 

Remark 4.6 We note that, to control the tail of the integral in (14.261 ) in d > 3 it suffices to 
invoke the diagonal heat-kernel estimate, 

E^(X t = x)<^, xeZ d , (4.29) 

which in the elliptic case is an immediate consequence of the mixing theory for Markov chains 
based on isoperimetric inequalities. This is sufficient because the finiteness of the Green function 
in d > 3 permits us to define (/, (— even for / > which enables us to reduce the 
general case to positive / by decomposing the test function into a positive and negative part and 
applying 

E v {f,^f)< H/ll?^, (4.30) 

which is uniformly integrable when d > 3. However, to include d = 2, we cannot disregard the 
cancellations due to the vanishing of f f(x)dx and thus the stronger derivative bound (13.231 ) is 
necessary. A similar situation occurred in Giacomin, Olla and Spohn ifTBl where a stronger Nash 
continuity estimate was required to include d = 2. 

Finally, we are ready to establish the main result of this paper: 

Proof of Theorem |2~?1 Let \i be a translation-invariant, ergodic gradient Gibbs measure with zero 
tilt and let jl be its extension to M B ( z<i ) x M B ( z<i ). We want to prove that 4>(f € ) tends weakly to a 
normal random variable with mean zero and variance 

(/, (-Qr 1 /)- By Lemma O it suffices 

to prove this for / G C$°(R d ) n H. 

By Lemma [3~4l we know that 4>(f ) is Gaussian conditional on k. The standard formula 

E{e x ) = e iE (*)-f v *-P0 (4.31) 
valid for any Gaussian random variable implies via ( I3.15I4XT61 ) 

E^e^l^Xn) = e -5(U-£«)- 1 f). (4.32) 
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By Proposition 14.41 we have (/ e , (— C K )~ l f e: ) —>■ (/, (— Q)~ x f) in //-probability. Since the 
right-hand side of (14.321) is a bounded continuous function of this inner product, (12.151 ) follows 
by means of the Bounded Convergence Theorem. □ 



5. Potential theory for random conductance model 

The proof of Lemma l331 leads us to the study of potential theory for operators depending on a ran- 
dom environment that belongs to the class of random conductance models. A good deal of what 
is to follow exists explicitly, or implicitly, in the literature. We borrowed some of the notation 
from the paper of Mathieu and Piatnitski [21] although the formalism draws on earlier works in 
homogenization theory; see, e.g., the book by Jikov, Kozlov and Oleinik iflTl . Notwithstanding, 
the content of Sect. I5.2l is presumably new. 

5.1 Basic notions. 

Consider a translation-invariant v probability measure on fi = M + satisfying the ellipticity 
condition 

3e>0: z/(e < K b < V e ) = 1, b G B(Z d ). (5.1) 
Let L 2 {u) denote the closure of the set of all local functions in the topology induced by the inner 
product 

{h, g) :=E V (h(K)g(K)). (5.2) 

Let B := {Si, . . . , e^} denote the set of coordinate vectors in 7L d . The translation by the vectors 
in B induce natural unitary maps Ti, . . . , on L 2 (v) defined via 

(Tjh) :=ho T ~ ej , j = 1, . . . , d. (5.3) 

Apart from the square integrable functions we will also need to work with vector fields by which 
we will generally mean measurable functions u: Q x B — > M.or Q x B —> R. d depending on the 
context. We will sometimes write u\, . . . , it^ for u(-,6i), . . . , u(-, e^) — note that these may still 
be vector valued. 

Remark 5.1 While we index vector fields only by the positive coordinate vectors, in certain 
situations it is convenient to have them defined also for the negative coordinate directions via 

u(k, -b) := -u(r^ b K, b), b G B. (5.4) 

As we will see this definition will automatically ensure that the cycle condition (see Lemma [5721 
below) holds for the trivial cycles crossing only a single edge. 

Let L 2 ec (u) be the set of all vector fields with (u,u) < oo, where (•, •) denotes the inner 
product 

(u, v) := E v ( K b u(k, b) ■ v(k, bU . (5.5) 
beB 

Examples of such functions are the gradients V/i of local functions h G L 2 (fi) defined compo- 
nentwise via the formula 

(V% := Tjh -h, j = l,...,d. (5.6) 
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We denote L^(u) the closure of the set of gradients of local functions in the topology induced by 
the above inner product. 

Lemma 5.2 Let u £ Ly. Then u satisfies the cycle condition 

it 

^ u{t X] h, Xj+i - Xj) = (5.7) 
3=0 

for any finite (nearest-neighbor) cycle {xq, x\, . . . , x n = xq) on Z, d . In particular, there exists a 
shift-covariant function u: Q x 7L d — > M. d such that u(k, b) = u(k, b)for every b € B. 



Proof. The cycle condition (15.71 ) holds trivially for all gradients of local functions. Indeed, if 
u = V/i then, in light of (15.41 ) we have 

u(T Xj K,Xj + i - Xj) = (Vh) Xj+1 -. Xj (T x .K) = ho t Xj+1 (k) - ho T Xj {n). (5.8) 

A corresponding limit extends this to all of Ly . To define u(-, x), we integrate properly shifted 
values of u along a path from zero to x; the cycle condition guarantees that the result is indepen- 
dent of the choice of the path and that u is shift covariant. □ 

We will henceforth use the convention to write u for the extension of a shift-covariant vector 
field u € Ly ec to a function on Z d . Notice that the shift Tj extends naturally via 



Tju(n,x) :=w(rg.«, x) = (Tju)(n,x). (5.9) 
Next let us characterize the functions in (Ly) : 

Lemma 5.3 For u E Ll ec (v), let Cu be the function in L 2 (v) defined by 

{Cu)(k) := ^2 [ K b u(K, b) - {r_ b K) b n(r_ b K, b)] , (5.10) 
beB 

where —b is the coordinate vector opposite to b. Then we have 

nG(L|) ± Cu = 0, u-a.s.. (5.11) 
Ifu satisfies the cycle condition and u is its extension, then Cu(t x k) = C k u(k, x). 

Proof. These are direct consequences of the definitions, the translation invariance of v and a 
simple calculation. □ 

Note that Cu plays the role of divergence — i.e., the total flow out of a given vertex — of vector 
field u. However, we prefer to denote it by C to emphasize its connection with the operator C K . 

5.2 Uniqueness of harmonic embedding. 

Clearly, all u G L v are shift covariant and have zero mean. A natural question is whether every 
shift-covariant zero-mean u is in Ly. (Note that this is analogous to asking whether every closed 
differential form is exact.) Our answer to this is affirmative: 

Theorem 5.4 Suppose v is ergodic. Then every u £ Ly ec (v) which obeys the cycle condition 
(15.71 ) and E v u = satisfies u € Ly. 
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Recall again that ( 15.71 ) and zero expectation are necessary for u 6 L v . The above says that 
these conditions are also sufficient. To prove the theorem, we will need the following lemma: 

Lemma 5.5 Let Pj denote the orthogonal projection onto Ker(l — Tj) in L 2 ec (v). Ifv is ergodic 
and u € L 2 e( .{y) satisfies (15.71 ). then PjU = E u (Pju), v-a.s.. 

Proof. Fix u € L 2 ec (u) that obeys ( 15.71 ) and let u be the corresponding shift-covariant function. 
We will prove the claim only for the component u\ = ei); the other cases follow analogously. 
By translation covariance and the L 2 -Ergodic Theorem, 

3^}Ell = iy T k Ul _^ P inL 2 (u). (5.12) 

k=0 

If v were separately ergodic (i.e., ergodic with respect to T\ alone) then the claim would imme- 
diately follow by the fact that every T% -invariant function must be constant. To make up for the 
potential lack of separate ergodicity, we note that translation covariance of u and the fact that u 
obeys the cycle conditions yield 

Tju(-, nei) = n(-, nei + e,-) — u(-, e,) 

- , - x , . J, , - x (5-13) 
= u{-,nei) - u(-,ej) +T nh u{-,ej) 

It follows that Tj^u(-,nei) also converges to P\U\ (in L 2 ) and so, by the continuity of Tj, 

T j Piu 1 = P 1 u 1 . (5.14) 

Hence P\u is invariant with respect to all shifts and so it is constant i^-a.s.. □ 

Proof of Theorem 15.41 Suppose v is ergodic and let u € L 2 ec (^) obey ( 15.71 ) and E v u = 0. 
The boundedness of the k&'s away from zero and infinity ensures that u € L 2 ec (z^) iff all of its 
components are in L 2 (u). Our task is to construct functions h e G L 2 (v) such that Vh t — > u 
in L 2 ec (i/). We define 

^ 1 + e n+1 

n>0 V ; 

and note that this is the unique solution of the equation (1 + e — T\)h e = — u\. This observation 
implies 

(1 - Tijfce = -ui - e/t e (5.16) 

and so the first component of V/i e converges to that of u provided eh e — > in L 2 (v). To see what 
happens with the other components of S/h e , we note that the cycle condition ( 15.71 ) translates into 

{1-T j )u 1 = (l-T 1 )u j . (5.17) 

Applying this to the definition of h e , we conclude 

(l-Tj)he = -Uj-eh e , (5.18) 

where h e is defined as h e but with u\ replaced by Uj. Again, it suffices to show that eh e — > 
in L 2 {y) which will boil down to the same argument as for j = 1. 
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To prove that eh e Owe note that, for the inner product in ( 15.21) , 

n + 1 



n>0 

Introducing the notation A n u for the average 



^ h *) = £ 7rT7^T2 (ni ' T ? Ui) - (5 - 19) 



1 " 

A n u := - Vlfu, 



n 

fc=0 



(5.20) 



inserting this into the above sum and reordering the terms, we get 

{he, he) = %^ + £ " 3 (e(n + 1) - 1) (tn^tn). (5.21) 

n>l ^ ' 

By Lemma 1531 the L 2 -Ergodic Theorem and the fact that u has zero expectation in v, 

A n u\ — ► P\U\ = E v u\ = 0, inL 2 (^), (5.22) 

n— too 

and so (ui, A n u\) — > as n — > oo. A straightforward estimate now shows that the sum in (15.211 ) 

is o(e~ 2 ) and so e 2 (h e , h e ) — > as desired. □ 

Not every function in L 2 ec (z^) necessarily belongs to Ly(i^). A prime example is the position 
vector field sc(«, 6) = b. Indeed, let \ '■ M B(z<i) xB^l d be the projection 

X := -proj L 2 (u) x. (5.23) 



Since % £ -^v(^) ^ satisfies (15.71 ) and we may extend it to a function K B ( Z ) x Z d — > R d 
by setting x(-,0) = an d integrating the gradients along oriented paths. Lemma 1531 implies 
that x+x is harmonic, C K (x+x) = 0. Moreover. Lemma l5.5l and (15.121 ) show that y(-. nej) /n — > 
and so x + x 7^ 0. It follows that x L v and so L v 7^ £yec- 

The function x is generally referred to as the corrector because it corrects for the non-harmoni- 
city of the position function. The corrector can be defined by appealing to spectral theory (Kipnis 
and Varadhan [18); see also Berger and Biskup [4]); the above "projection" definition is inspired 
by those in Giacomin, Olla and Spohn lPT3l or Mathieu and Piatnitski lf2TTl . As a side remark we 
note that the function x + x actually allows us to characterize the space of all square-integrable 
shift covariant functions: 

Corollary 5.6 Suppose v is ergodic. Then every shift-covariant M.-valued u G {L\) 1 ~ can be 
obtained from x+x, where x is the position function and — x is its orthogonal projection onto L v , 
by means of the linear transformation 

d 

Uj{K) = ^2[&k- {xj + Xj)]E u u k , j = l,...,d. (5.24) 

k=l 

Here u\, . . . ,Ud stand for u(-, 61), . . . , u(-, ej) and similarly for (M. d -valued objects) xi,...,Xk 
and Xli ■ ■ ■ j Xd- In particular, for "^-valued vector fields we have 

{u £ Ll ec (u) : shift covariant} = L% {A • (x + x) ■ A G R d }. (5.25) 
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Remark 5.7 The observation ( 15.251 ) — which in the language borrowed from differential ge- 
ometry implies that linear transforms of x + \ WQ the only closed forms that are not exact — 
has previously been discovered in the context of multicolored exclusion processes (Quastel ll24l 
Theorem 9.1]). However, while our proof is based only on soft, Poisson-equation based argument 
(15.15I45I221 ). that of [24] requires an explicit bound on the spectral gap for the corresponding dy- 
namics. The reason may be the reliance of j24l on the spatial ergodic theorem which, naturally, 
leads to bounds involving the Poincare inequality and/or spectral gap. 

Proof of Corollary 15.61 This is a simple consequence of part (1) of Theorem I54t Let w = 
{w\ , . . . , Wd) be the (R- valued) vector field defined by 

d 

wj := ^[efc • (xj +Xj)]E v u k , j = l,...,d. (5.26) 

k=l 

We will show that w = u. First, both u and w obey the cycle condition and so does u — w. As 
X G Ly{y) implies E v \ = the fact E v x k = x^ = &k shows E v Wj = E v Uj, i.e., E v (u — w) = 
0. Theorem l54[ l) implies u — w G Ly. On the other hand, C(x + x) = implies Cw = and 
so, by Lemma [531 w € (Ly)^. Thus u — w € (Ly) 1 - holds as well. It follows that u = w as 
claimed. 

As to (15.251 ). the argument we just used ensures that a shift-covariant u £ L^ ec (u) can be 
written as A • (x + x), where Xj := E u Uj, plus a shift-covariant vector field with zero expectation. 
By Theorem EH 1 ), the latter is in Ly. □ 

We still owe the reader the proof of Lemma [331 

Proof of Lemma [PI Since g has square-integrable gradients, we have g G L^ ec (z^). As g is shift 
covariant and has zero expectation, Theorem [53] implies that g € Ly. But g is also harmonic and 
so Lemma [531 in turn forces g G (Ly) . Thus g = as desired. □ 

Recall that % is the extension of x subject to the condition 0) = 0. With this object at 
hand, we may even remove the restriction to zero slope in Lemma [34] 

Corollary 5.8 Let fibe a translation-invariant, ergodic gradient Gibbs measure and let jl be its 
extension to x R B ( zd ). Let t G R d be the tilt of fi and let & := cr({K b : b G B(Z d )}). 

Then /2(— \JP)(k) is Gaussian with mean 

Efi(fa-<fo\&)(K)=t' [x + x{k,x)] (5.27) 

and the covariance given by (— £ K ) _1 . 

Proof. Let u denote the extension of the shift-covariant function £^ (770,6 b G B. Clearly, 
we have 

u(k,x) =£/i(0x- 00 ooeZ d . (5.28) 
Inspecting the proof of Lemma [341 all formulas carry over except (13.221 ) which becomes 

Efi(u{-,x)) =t- x, x£Z d . (5.29) 

Thus u — t ■ (x + x) is harmonic, shift-covariant, and of zero mean and so u = t ■ (x + x) by 
Lemma [531 and Theorem 1541 □ 
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6. Discussion and open problems 

The proofs in the present note often rely on the fact that the random walk with generator C K is 
uniformly elliptic. This enters via the assumption that g, defining the potential V, has support 
bounded away from zero. While we believe that the general picture carries over even if we let 
the support extend all the way to zero, a number of steps in the proof become quite subtle. For 
instance, the pointwise heat-kernel asymptotic for this walk may take a radically different form 
(Berger, Biskup, Hoffman and Kozma [3]) and it is not known under what conditions on the en- 
vironment the walk scales to Brownian motion. Progress in this direction for i.i.d. (or i.i.d. domi- 
nated) environments has been made only recently (Mathieu [20], Biskup and Prescott [6], Barlow 
and Deuschel Q). 

Another interesting open question concerns the annoying restriction to zero-tilt gradient Gibbs 
measures. As is seen from Corollary 15.81 once the tilt is nonzero, (14.321 ) has to be modified into 

£_( e #(/)-% W))] | = e -i (W)-f (/,(-£«)- x /). (61) 

When we plug in f e for /, the second term in the exponent still converges to (/, (— Q) -1 /) in 
probability as e J. 0. To show convergence of cj) e — E^(\> e to Gaussian free field, we thus have 
to show that [t ■ x)(fe) converges in law to a normal random variable. If V is strictly convex 
and of the form (11.41 ). we in fact know this to be true by the main result of Giacomin, Olla and 
Spohn [13]. However, the general case seems to be a hard open problem (cf. [4, Conjecture 5]). 
Already in d = 1, where the corrector is simply the sum 

i x ~ 1 i 

= pi £( C), (6.2) 

with C := E^l/ttft), the desired result requires that the resistivities 1/kj, exhibit a specific 
correlation decay ifTTl Theorem 7.6]. Very little information concerning this is known in d > 2; 
see Mourrat 11221 for some recent progress in this direction. 

Our final remark concerns a generalization to potentials V where (11.41) has been modified into 

V(rj) := - log J Q{dHc)e~ w ^\ (6.3) 

with (W K ) denoting a family of strictly convex, even, measurable functions with uniformly su- 
perlinear growth at ±oo and a uniform lower bound. In this case we may still consider the ex- 
tended gradient Gibbs measures; however, the conditional law given the re's is no longer Gaussian. 
Notwithstanding, given k, the Helffer-Sjostrand representation still applies and leads to a random 
walk in a dynamic random environment. Its annealed CLT would imply the Gaussian scaling 
limit for the 0-field. 
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